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摘 要 : 卫星 对 地 观测 及 导航 定位 的 快速 发 展 对 地 球 同步 轨道 (Geostationary Earth Orbit, GEO) 
卫星 的 定 轨 精度 有 了 米 级 其 至 更 高 的 需求 。 本 文 研究 利用 其 长 基线 干涉 (Very Long Baseline 
Interferometry, VLBI) 时 延 数 据 进行 GEO 卫星 快速 定位 的 方法 。 首 先 对 VLBI 时 延 进行 理论 精度 分 析 ， 然 后 

利用 中 国 VLBI 网 进行 GEO 卫星 观测 试验 及 数据 相关 处 理 ， 最 后 利用 短 时 VLBI 时 延 数据 完成 对 GEO 
卫星 的 定位 计算 。 针 对 时 延 中 存在 与 测 站 和 基线 相关 的 异 党 常 值 问题 ， 本 文 提 出 了 相应 的 检测 校正 方法 。 数 
据 经 校正 后 重 新 处 理 ， 得 VLBI 残余 时 延 测量 精度 约 0.16ns， 等 效 于 约 0.58 米 的 卫星 横向 位 置 误差 。 
2 I GEO 卫星 的 快速 高 精度 定位 及 轨道 快速 恢复 。 
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地 球 同 步 轨 道 (Geostationary Earth Orbit, GEO) 卫星 轨道 高 约 36000 公里 ， 绕 地 公转 周期 与 地 
球 自转 周期 一 致 ， 一 定 范围 内 相对 地 面 点 静止 ， 在 大 地 测量 、 导 航 定位 和 灾害 监测 等 方面 发 挥 重要 作 
用 。 由 于 GEO BE 运动 中 受 各 种 摄 动力 影响 ， 特别 是 其 太阳 能 避 板 受到 的 太阳 光 压 摄 动 大 ， 使 卫星 
的 实际 轨道 发 生 漂移 ， 偏 离 平 衡 点 中。 为 减少 轨道 误差 对 用 户 定 位 精度 的 影响 ， 需 尽快 获取 卫星 准确 
位 置信 息 ， 因 此 研究 GEO 卫星 机 动 后 的 位 置 快速 测量 方法 具有 重要 意义 。 
目前 对 ba 卫星 视 向 位 置 测 量 主要 依靠 测 距 技术 ， 包 括 伪 距 和 载波 相位 测量 、 转 发 式 测 距 测量 
和 激光 测 距 等 。 高 精度 测 距 需 卫星 转发 测 距 信 号 ， 为 主动 式 测量 ， 定 轨 精 度 百 米 量 级 外 。 制 约 测 距 定 
轨 精 度 的 因素 与 GEO 卫星 特点 有 关 : 卫星 轨道 高 ， 地 面 跟踪 站 构 网 范围 相对 较 小 ， 导 致 观测 网 几何 
结构 略 差 ， 相 对 静止 特性 使 得 星 站 间 动 力学 和 几何 约束 变 弱 ， 造 成 钟 差 和 测 站 偏差 等 系统 误差 解 算 困 
XE; 为 维持 同步 轨道 ， 卫 星 需 频繁 轨 控 ,但 星 上 发 动机 推力 不 易 建 模 ， 造 成 后 续 轨 道 恢复 和 预测 困难 


对 GEO 卫星 轨道 平面 位 置 测 量 可 采用 光学 成 像 或 射电 干涉 等 被 动 式 测量 方法 ， 这 些 方法 优点 为 
仅 需 接收 卫星 下 行 信 息 ， 不 占用 星 上 资源 。 拥 有 独立 本 振 的 其 长 基线 干涉 (Very Long Baseline 
Interferometry, VLBI) 技术 是 一 种 射电 干涉 测量 技术 ， 因 其 具有 超 高 测 角 分 辨 紊 和 灵敏 度 ， 已 被 广 
泛 应 用 于 深 空 探 测 器 跟踪 测量 门 。VLBI 对 确定 卫星 横向 (垂直 于 卫星 轨道 径 向 平面 ) 位置 具 有 高 灵 
敏 度 ， 可 提供 观测 目标 角 位 置 的 几何 约束 。 差 分 VLBI 技术 通过 快速 交 蔡 观测 目标 探测 器 和 邻近 致密 
射电 源 ， 利用 信号 差分 可 有 效 减 少 介 质 时 延 、 仪 器 时 延 、 测 站 钟 差 等 测量 误差 ， 测 角 精 度 优 于 毫 角 秒 

联合 使 用 VLBI 技术 与 测 距 技 术 可 进一步 提高 GEO 卫星 的 三 维 位 置 解 算 精 度 。1984 年 Tadashi 
Shiomi 采用 差分 VLBI 技术 与 无 线 电 测 距 对 GEO 卫星 进行 联合 定 轨 试验 得 到 卫星 轨道 精度 约 100 米 
中 。 近 年 来 国内 学 者 也 开展 了 对 GEO. 卫星 的 仿真 试验 和 实测 研究 。2003 年 舒 着 春 利用 双 差 单 向 测 距 
(Delta Differential One-way Ranging, ADOR) 和 双 差 单 向 测速 两 种 差分 VLBI 技术 ， 开 展 我 国 首 次 
GEO 卫星 的 VLBI 实测 试验 ， 获 得 卫星 信号 干涉 条 纹 ， 未 进行 实际 卫星 定位 定 轨 计 算 。 试 验 和 DOR 
总 误差 约 41cm， 理 论 上 相当 于 用 2000km 基线 观测 地 球 同 步 轨道 目标 时 gm 的 位 置 误差 路。2006 年 杜 
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兰 进 行 数值 仿真 试验 表明 ， 
未 开展 实测 试验 。 
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在 测 距 观测 中 增加 VLBI $ 


2011 年 黄 勇 利 月 


BEDS. SEE 


H C 波段 转发 式 测 


LARK GEO 卫星 轨道 精 
距 数 据 和 VLBI 时 延 数 据 进行 


ChinaXiv 合 ; 


综合 定 轨 试 验 ， 定 轨 后 时 延 数 据 残 差 均 方 根 值 约 3.6ns， 时 延 率 约 0.4ps/s， 相 当 于 轨道 横向 位 置 误差 


约 10mU?!, 


2017 *E 5 SIZE RH 


0.14ns 的 时 延 


鉴于 
究 GEO 


卫星 进行 定位 计算 。 


EGRE, [n] 
前 述 试验 结果 均 未 与 观测 目 
卫星 轨 控 后 轨道 恢复 时 间 长 、 
卫星 位 置 的 快速 高 精度 测量 技术 ， 具 
VLBI network, CVN) 所 属 4 个 台 站 观测 GEO 卫星 ， 并 利 月 
CVN 由 天 马 (Tm) 、 


标 事 后 精 轨 比较 ， 淮 


北京 Bj) 、 


重要 现实 意义 。 我 们 组 乡 


昆明 (Km), 5 


时 间 段 重复 定 轨 结 果 的 中 误差 在 米 量 级 050， 为 结 
E 确 度 暂 无 法 衡量 。 
WEK, HX GEO 卫星 定 轨 精度 的 要 求 越 来 越 高 ， 


日 宽带 VLBI 系统 与 轨道 跟踪 和 转移 测定 系统 观测 GEO 卫星 ， 得 到 约 


果 的 内 符合 精度 评估 。 
因此 ， 研 


! 中 国 VLBI 网 (Chinese 
H VLBI 短 时 观测 的 时 延 测 量 数据 对 GEO 
鲁 木 齐 (UD 4 观测 站 和 位 


于 上 海 的 数据 处 理 中 心 组 成 ， 含 6 条 观测 基线 。 本 文 首先 从 理论 上 对 GEO 卫星 VLBI 时 延 观测 精度 


进行 分 析 ， 然 后 利用 约 4 小 时 的 S 波段 VLBI 时 延 数 据 进行 卫星 定位 计算 。 同 时 文中 详 


描述 了 


VLBI 时 延 结果 中 的 异常 值 检测 和 校正 方法 ， 并 利用 校正 后 的 VLBI 时 延 数据 和 新 开发 的 后 处 理 程序 


重新 进行 GEO 卫星 定位 计算 ， 


通过 定位 结果 验证 时 延 跳 变 的 校 了 


E 效 果 。 将 计算 结果 与 国际 全 球 导航 


卫星 系统 服务 组 织 CInternational Global Navigation Satellite Systems Service, IGS) 提供 的 事后 精 轨 数 
据 进 行 比较 ， 可 作为 定位 结果 的 外 符合 准确 度 估计 。 


1 卫星 差分 VLBI 观 测 原 理 与 精度 分 析 


1.1 


的 角 距 范围 越 入 。 
。 图 中 1、 


通过 两 
示 2 为 两 观测 台 站 ，B 为 两 台 
台 站 观测 卫星 方向 矢量 ，7 为 射电 源 观 测 时 延 。 


卫星 差分 VLBI 观测 原理 


VLBI 观测 量 为 时 延 和 时 延 率 ， 即 两 个 地 面 天 线 接收 同一 
括 群 时 延 〈 率 ) 和 相 时 延 CK) 。 
作为 参考 源 ， 交 蔡 观 测 参考 源 及 目标 探测 
言 号 差分 ， 可 有 效 减 少时 延 测 
站 组 成 的 基线 矢量 , ?为 台 
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差分 VLBI 技术 即 选 取 目 


音 号 的 时 间 延 迟 和 时 间 延 迟 变化 率 ， 包 


m E o 射电 源 


标 探 测 器 附近 位 置 精 
器 。 一 般 要 求 两 者 角 距 小 于 10° 0 
量 误差 。 差分 VLBI E 图 如 图 


可 用 
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站 观测 射 
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1 差分 VLBI 原理 示意 区 
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Fig. 1 Schematic 


iagram of Differential VLBI Principle 


测量 时 ， 按 照 “参考 源 -目标 卫星 -参考 源 ” 的 顺序 观测 ， 


测量 | 5&4 (t2). 在 t。 时 刻 ， 


Eti» t2» 
为 T9(t1)，T5(t2)，T94(t3)， 上 标 q 表 示 参 考 源 观测 ，s 表 示 卫 星 观 测 。 通 
有 如 下 关系 


T(t2) = TY (t2) + Th (tz) 


T(t2) = Tg(t2) 十 Terr(tz) 


2 


ts 时 刻 对 应 的 VLBI 时 延 值 分 别 
过 插值 可 得 到 ts 时 刻 参 考 源 观 


(1) 


(2) 


(FEAF 


) OOO 


其 中 ，rg 为 几何 时 延 ， 可 以 由 台 站 坐标 、 观 测 源 坐标 和 地 球 定 向 参数 等 先 验 值 计 算得 到 。rerr 为 
观测 中 的 系统 误差 和 随机 误差 项 之 和 ， 包 括 介 质 时 延 、 钟 差 、 仪 器 时 延 和 观测 噪声 等 。 差 分 VLBI 时 
延 观测 值 为 


AT = T(t2)— T1 (t2) = 15(t2) — Th (t2) 十 Er (3) 


Er = Têrr(t2) — Tarr (t2) (4) 


Alti (t2) c2, (t2)23 8] — E Zap d vx ze de PAARE ESEE, A E A 
时 间 短 且 参 考 源 和 目标 探测 器 角 距 近 ， 所 以 介质 时 延 相差 小 。 因 此 ， 差 分 得 到 的 ei 为 极 小 误差 项 ， 消 
除了 VLBI 观测 中 无 法 被 精确 建 模 的 共同 误差 源 影响 ， 得 到 卫星 的 实际 VLBI 时 延 值 Tr5(t,) 


t*(t2) ~ 1*2 (t5) + AT (5) 
zs(tz) 包 括 的 卫星 的 位 置信 息 ， 可 用 于 定位 。 
1.22 ”时 延 精度 与 卫星 位 置 误差 分 析 


为 评估 卫星 位 置 解 算 精 度 ， 需 对 VLBI 测量 精度 进行 理论 分 析 。VLBI 时 延 测 量 误差 主要 包括 射 
电源 源 位 置 、 台 站 位 置 、 地 球 自转 参数 和 介质 等 带 来 的 系统 性 误差 ， 以 及 观测 源 噪 声 、 氢 原子 钟 不 稳 


定 度 、 介 质 扰动 和 仪器 相位 拌 动 等 带 来 的 随机 误差 "”。 
射电 源源 位 置 误差 VLBI 观 测 时 从 国际 天 球 参考 架 优先 选择 位 置 精度 高 的 点 源 ， 因 此 系统 误差 
中 的 源 位 置 时 延误 差 是 一 个 与 源 本 身 位 置 精度 和 投影 基线 长 有 关 的 极 小 量 。 估 算 公 式 为 上 


Em = E (6) 


B, 为 基线 在 天 空 平面 投影 的 长 度 ，eg 为 河 外 源 位 置 误 差 ，c 为 光速 。 大 部 分 源 的 位 置 精 度 在 0.5- 
2.0nrad， 对 投影 长 度 3000km 的 基线 ， 射 电源 位 置 误差 0.75nrad 时 ， 由 射电 源 位 置 误差 带 来 的 时 延误 
差 约 为 0.0075ns。 

台 站 位 置 误差 ” 台 站 位 置 时 延误 差 可 用 基线 在 天 球面 的 投影 误差 表示 ， 与 射电 源 和 卫星 夹 角 有 
关 ， 表 示 为 上 


EAT F = (A6) €BL (7) 


A0 735; FA URURIL D EE SER, eg, 为 基线 投影 误差 ， 假 定 射 电源 和 卫星 夹 角 0.2rad， 基 线 投影 误差 
0.02 米 ， 则 由 台 站 位 置 误差 带 来 的 时 延误 差 约 为 00067ns。 

地 球 定向 参数 误差 地球 定向 参数 对 时 延 的 影响 主要 由 模型 计算 时 坐标 系 转换 时 引入 ， 取 决 于 
UTI 和 极 移 参数 的 精度 ， 与 卫星 和 参考 源 间 夹 角 成 比例 ， 关 系 为 上 9 


1 
EAr F E (A0) surpw (8) 


surpw 为 地 球 定向 参数 UTI 和 极 移 的 误差 ， 为 一 极 小 值 ， 假 定 射电 源 和 卫星 夹 角 0.2rad， 球 定向 
参数 为 0.02 米 ， 则 由 其 产生 的 时 延误 差 约 为 0.0067ns。 

中 性 大 气 系统 误差 ” 该 误差 与 台 站 天 顶 方向 的 干 湿 对 流 层 误差 及 台 站 观测 目标 高 度 角 有 关 ， 近 
似 关 系 如 下 


pz 1 1 
Eaa er (9) 
c |sinys + 0.015 siny; + 0.015 


分 别 为 台 站 观测 卫星 和 射电 源 的 高 度 角 ，p;z 为 天 顶 方向 的 干 湿 对 流 层 误差 。 假 设 卫星 高 度 

fg 40”， 射 电源 高 度 角 39”， 天 顶 方向 的 干 对 流 层 误差 为 0.002m， 天 顶 方向 的 湿 对 流 层 误差 为 
0.005m， 两 个 台 站 干 湿 对 流 层 影 响 产 生 的 时 延误 差 约 为 0.0075ns。 
电离 层 系 统 误 差 ”电离 层 系 统 性 误差 对 时 延误 差 的 影响 较 复杂 ， 可 由 Klobuchar 模型 简略 计算 ， 

若 观测 台 站 周围 布设 有 双 频 GPS 接收 机 可 获得 更 准确 的 台 站 天 空 电 离 层 信 息 。VLBI 观测 电离 层 时 延 
误差 影响 的 典型 值 为 0.02ns。 
卫星 热 噪声 误差 GEO 卫星 热 噪 声 与 卫星 信号 信 噪 比 、 观 测 通道 带宽 和 积分 时 间 等 有 关 ， 近 似 
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关系 为 


VE 
Pr 2nAfow,SNRs 


Afsw 为 卫星 观测 带宽 ，SNRs 为 卫星 信号 信 品 比 。GEO 卫星 数 传 信号 通常 为 带宽 10MHz 的 单 通 
道观 测 ， 假 定 卫 星 信 噪 比 1000， 则 卫星 热 噪声 带 来 的 时 延误 差 约 0.0551ns。 

台 站 钟 不 稳定 度 误差 ” 氧 原子 钟 不 稳定 性 带 来 的 误差 具体 可 用 目标 交替 观测 间隔 内 的 台 站 钟 频 
率 稳定 性 表示 ， 公 式 为 


(10) 


En = Ts_qéEAf/f (11) 


Ts_o 为 射电 源 和 卫星 交 蔡 观测 时 间 ，enypjr 为 站 钟 的 稳定 性 ， 若 探测 器 和 河 外 源 观测 时 间 间 隔 为 
600s， 钟 频率 稳定 度 为 10-14， 则 其 时 延 影响 为 0.0060ns。 
仪器 相位 抖动 误差 ”仪器 相位 抖动 产生 的 时 延误 差 与 信号 总 带宽 近似 成 反比 ， 关 系 为 0 


Enr = V2V2— (12) 


mT 


so 为 仪器 相位 误差 ， few, 为 总 观测 带宽 。 假设 仪器 相位 抖动 为 0.2”， 总 观测 带宽 为 40MHz， 则 
时 延误 差 为 0.0278ns。 
介质 扰动 误差 对流层 和 电离 层 的 随机 性 误差 与 各 自 观测 方向 的 扰动 不 确定 性 及 卫星 和 参考 源 
间 夹 角 有 关 ， 该 项 不 能 通过 差分 运算 消除 ， 表 示 如 下 中 
1 A8 
EAT F z 00745) Pfluct (13) 


prmct 为 10” 分 离 角 的 对 流 层 或 电离 层 波动 不 确定 性 ， 若 对 流 层 和 电离 层 波动 不 确定 性 均 为 
0.01m， 则 两 者 对 时 延 影响 均 为 0.0191ns。 

上 述 各 误差 项 的 统计 结果 见 图 2 ， 图 中 红色 项 为 系统 误差 ， 蓝 色 项 为 随机 误差 ， 黑色 项 为 所 有 误 
差 项 平方 和 的 平方 根 值 (Root Sum Square, RSS) ，RSS 总 误差 约 0.0720ns。 
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图 2 时 延误 差 估计 


Fig. 2 Estimation of Delay Error 


差分 VLBI 时 延 测量 误差 er 与 卫星 垂直 于 径 向 轨道 的 平面 位 置 误 差 5 之 间 满 足 " 
ôs B 
€ 


EAr F 


:一 ( 14) 
L 


式 中 OB 为 基线 在 天 空 平面 投影 的 长 ,区 近似 认为 是 测 站 到 卫星 距离 。 对 于 GEO 卫星 , LAN 
36000km， 若 投影 基线 长 3000km， 差 分 VLBI 时 延 测 量 精 度 在 0.0720ns 时 ， 轨 道 横向 位 置 精 度 约 
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0.26 米 的 。 

为 了 对 定位 结果 的 可 靠 性 进行 直观 评估 ， 
可 将 试验 结果 与 IGS 提供 的 事后 精 轨 数 据 进行 
比较 。 不 同 机 构 提 供 的 GEO 精 轨 结果 本 身 存 
在 一 定 差异 。 目 前 ， 德 国 地 球 科学 研究 中 心 
CGeoForschungsZentrum, GFZ) 和 武汉 大 学 
(Wuhan University, WHU) 两 家 机 构 可 公开 
提供 观测 日 滞后 时 间 约 十 几 天 的 GEO 卫星 精 
轨 。 两 机 构 提供 的 观测 日 北斗 导航 C03 星 精 轨 
比较 结果 如 图 3 所 示 。 图 中 横 坐 标 表 示 观 测 日 
时 刻 ， 纵 化 标 表示 两 机 构 精 轨 坐 标 差 值 。 赣 楼 
绿 黑 散 点 分 别 为 X, Y, Z 坐标 分 量 差 和 三 维 
坐标 差 ， 均 方 根 值 依次 为 1.09m，0.64m ， 


这 


0.17m 和 1.27m， 三 维 坐标 差 最 大 值 约 为 1.45m。 


说 明 不 同 机 构 提供 的 精 轨 数 据 本 身 包含 米 量 级 
的 偏差 。 


2 实测 数据 处 理 与 分 析 


2.1 观测 设置 


我 们 在 UTC 时 间 2018 年 1 月 13 日 组 织 了 
在 2.217GHz， 具 体 观 测 参 数 设 置 见 表 1。 


表 1 观测 参数 设置 
Tab. 1 Setting of Observation Parameters 


code b8113a 


3:30 ~ 7:30 


1741-038, 1908-201, 2008-159 
2134400, 2227-088, 3C454.3 


Beidou C03 


Beijing (Bj) , Kunming (Km) 
Tianma (Tm) , Urumqi (Ur) 


S (2217GHz) 


observation period 
quasar targets 
satellite target 
stations 


beacon frequency 


channels 8 
bandwidth 4MHz 
sampling rate 8MS/sec 
quantization bits 2bit 
data recording format Mark5B 


CVN 的 4 台 站 参与 观测 ， 最 长 基线 为 东西 向 上 
最 短 基 线 为 北京 -天 马 基 线 (BjTmO ， 长 度 约 11 


coordinate ditterence (m) 
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图 3 GFZ 和 WHU 的 GEO 卫星 精 轨 比较 
Fig. 3 Comparison of Precise GEO Satellite Orbi 
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次 约 4 小 时 的 GEO 卫星 差分 VLBI 观测 ， 信 标 频率 
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图 4 射电 源 和 卫星 天 球 位 置 及 轨迹 
Fig. 4 Position and Trajectory of Radio Sources and 
Satellites in Celestial Reference System 


的 乌鲁木齐 -天 马 基 线 (Ur-Tm) ， 长 度 约 3249km; 
14km。 由 于 GEO 卫星 轨道 高 度 较 河 外 射电 源 低 ， 


且 相 对 射电 源 每 小 时 沿 天 赤道 方向 移动 15 度 左右 ， 导 致 不 同 测 站 观测 


星 和 邻近 参考 源 的 视差 角 不 


同 ， 且 大 小 随时 间 变 化 。 为 提高 观测 精度 ， 从 国际 天 球 参考 架 中 选择 6 颗 邻 近 C03 卫星 的 致密 射电 


源 作为 参考 源 ， 射 电源 和 卫星 天 球 位 置 及 轨迹 如 图 


2.2 数据 处 理 


4 所 示 。 


利用 CVN 软件 相关 处 理 机 scorr 进行 原始 观测 数据 的 相关 处 理 后 ， 经 校正 可 得 到 观测 时 延 [。 
相关 处 理 时 利用 预先 计算 的 模型 时 延 对 两 台 站 接收 的 两 路 信号 进行 补偿 ， 得 到 残余 时 延 拟 。 射 电源 
模型 时 延 为 远 场 模型 ， 卫 星 模型 时 延 为 近 场 模型 ， 卫 星 模型 时 延 根据 卫星 的 预报 轨道 获得 。 均 由 上 海 
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天 文 台 VLBI 中 心计 算 。 

相关 处 理 时 侍 里 叶 点 数 设 为 1024， 卫 星 频 域 512 个 数据 点 的 互相 关 功 率 相 位 谱 结果 如 图 5 所 示 。 
6 个 子 图 分 别 为 北京 -昆明 基线 (Bj-Km) 、 北 京 -乌鲁木齐 基线 (BUD 、 北 京 - 天 杞 基线 (Bj-Tm) ~ 
昆明 -乌鲁木齐 基线 (Km-Ur) 、 昆 明 -天 马 基线 (Km-Tm) 和 乌鲁木齐 -天 马 基线 (Ur-Tm) 结果 。 图 
中 横 坐 标 为 频率 ， 单 位 kHz; 纵 坐 标 蓝 色 为 幅度 谱 ， 单 位 dB, EENM, BANE COO. 
相关 处 理 结果 可 以 看 出 ，6 条 基线 均 获得 较 好 的 干涉 条 纹 。 
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图 5 6 条 基线 互相 关 功 率 相位 谱 


Fig. 5 Cross-correlation Power Phase Spectrum of 6 Baselines 


对 相关 处 理 后 时 延 数 据 进 行 拟 合 、 钟 差 改 正 及 介质 改正 等 后 处 理 的 结果 如 图 6 所 示 ， 横 坐标 为 观 
测 时 刻 ， 纵 坐标 为 6 条 基线 的 残余 时 延 值 ， 单 位 ns. 6 条 基线 残余 时 延 拟 后 残 差 中 误差 分 别 为 2.35ns， 
1.15ns，0.39ns，0.75ns，0.70ns，0.86ns， 均 值 为 1.03ns。 图 上 显示 部 分 时 间 段 的 残余 时 延 结 果 存 在 
缺 值 和 跳 变 情况 ， 与 台 站 数据 和 数据 处 理 中 引入 的 误差 有 关 。 残 余 时 延 结果 与 模型 时 延 的 和 为 总 时 延 
数据 ， 因 VLBI 对 确定 卫星 径 向 位 置 不 敏感 ， 若 仅 利 用 总 时 延 数 据 解 算 卫 星 位 置 ， 解 算 结果 精度 较 差 ， 
因此 可 采用 距离 约束 的 方法 。 假 定 卫 星 的 地 心 距 等 于 由 其 预报 轨道 计算 的 地 心 距 ， 则 地 心 距 约束 的 先 
验 误差 由 预报 轨道 误差 给 定 04。 将 该 等 式 与 VLBI 观测 等 式 共同 组 成 误差 方程 解 算 卫 星 位 置 ， 并 与 
IGS 提供 的 精密 轨道 数据 进行 比较 ， 三 维 坐标 差 结果 如 图 7 所 示 。 两 者 坐标 差 小 于 35 米 ， 坐 标 均 方 
根 误差 约 16.76 米 。 图 7 中 部 分 时 间 段 卫星 定位 结果 存在 跳 变 ， 与 残余 时 延 结 果 图 中 存在 跳 变 的 时 间 
段 较 一 致 ， 对 应 时 间 段 数据 分 别 在 两 图 中 由 蓝 红 绿色 圈 标 出 。 
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Fig. 6 Residual Delay of 6 Baselines Fig. 7 Satellite Coordinate Difference 
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23 ”残余 时 延 改进 方法 


时 延 异 常 值 将 直接 影响 卫星 最 终 定位 结果 。 为 提高 卫星 位 置 解 算 精度 ， 需 对 VLBI RRN EX 
异常 值 进行 检测 和 校正 ， 改 正 VLBI 残 余 时 延 跳 变 问题 。 时 延 异 常 值 产生 原因 主要 分 两 类 : Qf 
相关 ， 如 台 站 测量 记录 设备 引入 的 误差 或 钟 差 修正 残余 部 分 ，@ 与 基线 相关 ， 由 数据 处 理 误差 引入 。 
因 残 余 时 延 数 据 为 时 间 序 列 数据 ， 异 常 检测 时 不 能 作为 孤立 的 样本 点 处 理 ， 借 助 多 项 式 拟 合 方法 和 基 
线 理论 时 延 闭合 差 为 零 的 约束 条 件 来 检测 两 类 异常 值 。 将 每 条 基线 的 残余 时 延 拟 合 值 残 差 作为 异常 值 
检测 数据 ， 偏 离 平 均值 一 定 倍数 标准 差 的 数据 即 为 异常 值 。 与 台 站 和 基线 相关 的 异常 值 检测 校正 方法 
流程 图 分 别 见 图 8 和 图 9， 主要 步骤 包括 中 确定 标记 异常 值 ， 凶 判断 是 否 台 站 、 基 线 相 关 ;， 名 计算 校 
正 值 并 校正 。 


esiduals were calculated 
after fitting 


End of scan cycle 


NO 
Number of combinations in 


which the closure error exceeds 
the threshold 


Determine the baseline to be 
corrected 


Identify and flag outliers Calculate the baseline 


correction value and correct it 


图 8 台 站 异常 值 检测 校正 流程 图 9 基线 异常 值 检测 校正 流程 
Fig. 8 Outlier Detection and Correction Flow of Stations Fig. 9 Outlier Detection and Correction Flow of Baselines 


首先 检测 台 站 相关 异常 值 ， 试 验 中 采用 拟 合 值 残 差 1 倍 标准 差 作 为 异常 值 检测 阔 值 ， 若 同一 时 间 
段 与 某 一 台 站 相关 的 所 有 基线 均 标 记 有 异常 值 ， 则 认为 该 时 间 段 存在 台 站 相关 的 数据 异常 。 在 确定 为 
台 站 相关 异常 值 后 ， 对 异常 时 间 段 的 台 站 数据 进行 校正 : 在 0 到 最 大 异常 偏离 值 范围 内 ， 每 隔 一 较 小 
间隔 如 0.0lns 确定 一 改正 值 ， 计 算 改 正 值 改正 后 所 有 基线 的 时 延 测 量 值 与 拟 合 值 差 值 平 方 和 ， 最 终 
的 改正 值 应 使 该 平方 和 最 小 。 

针对 基线 的 数据 处 理 引 入 的 误差 ， 误 差 检测 与 校正 方法 主要 基于 多 基线 理论 时 延 闭 合 差 为 零 的 条 
件 ，4 个 台 站 6 条 基线 共 4 种 组 合 方式 ， 如 用 BK-KU-BU 表示 Bj-Km、Km-Ur 和 Bj-Ur 基线 组 合 。 计 
算 不 同 基线 组 合 的 时 延 闭合 差 ， 求 每 个 scan 内 每 一 时 延 闭 合 差 的 平均 值 (以 下 称 scan 闭合 差 均 值 ) 。 
BEN 1.25ns, scan 闭合 差 均值 超 冰 值 处 为 需要 校正 的 时 间 段 ， 取 scan 闭合 差 均 值 绝对 值 的 基线 
组 合 闭 合 差 之 平均 值 为 校正 值 。 正 负 号 据 scan 闭合 差 均值 正 负 和 该 基线 在 组 合 中 位 置 判 新 ， 例 如 : 
在 某 一 scan 内 ，BK-KU-BU 和 BK-KT-BT 基线 组 合同 时 被 标记 ， 则 共同 基线 BK (BH Bj-Km 基线 ) 
需要 被 校正 ， 若 该 scan 内 BK-KU-BU 和 BK-KT-BT 基线 组 合 闭 合 差 均值 都 为 正 ， 则 BK 基线 需 减 去 
改正 值 ， 反 之 需 加 上 改正 值 。 台 站 记录 数据 缺失 会 导致 基线 数据 空缺 问题 ， 受 部 分 时 间 段 数据 空缺 和 
不 同 基 线 组 合 闭 合 差 差别 较 大 等 的 影响 ， 可 能 导致 无 法 确定 需 校正 的 具体 基线 ， 可 通过 修改 闭合 差 立 
值 的 大 小 来 综合 判断 。 

校正 前 后 的 基线 组 合 时 延 闭 合 差 对 比如 图 10 所 示 ， 横 坐标 为 观测 时 刻 ， 纵 坐标 为 基线 组 合 的 时 
延 闭 合 差 ， 单 位 ns。 蓝 色散 点 为 校正 前 的 时 延 闭合 差 ， 部 分 时 间 段 闭合 差 结 果 较 差 ， 红 色散 点 为 经 
台 站 和 基线 异常 值 校正 后 的 闭合 差 结果 ， 校 正 后 的 结果 更 接近 理论 零 值 。 校 正 前 后 的 残余 时 延 分 别 为 
图 11 中 蓝 色 星 号 和 绿色 圆圈 ， 从 图 中 对 比 看 出 ， 残 余 时 延 值 经 过 校正 后 变 平滑 。 校 正 后 的 6 条 基线 
残余 时 延 拟 后 残 差 中 误差 分 别 为 0.56ns，0.53ns，0.31ns，0.66ns，0.50ns，0.70ns， 均 值 为 0.54ns。 
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为 进一步 提高 卫星 定位 精度 ， 对 数据 重新 处 理 后 ， 残 余 时 延 数据 变 得 更 加 平滑 ， 见 图 11 中 的 红 
色 线 。 采 用 五 阶 多 项 式 拟 合 VLBI 残余 时 延 数 据 ， 时 延 精度 用 拟 后 残 差 中 误差 表示 。 最 终 得 到 6 条 基 
线 的 时 延 拟 后 残 差 中 误差 分 别 为 0.16ns，0.23ns，0.14ns，0.14ns，0.14ns，0.13ns， 均 值 为 0.16ns， 等 
效 于 轨道 横向 位 置 误差 约 0.58 米 。 

分 别 用 校正 后 和 重新 处 理 的 总 时 延 数据 对 GEO 卫星 定位 计算 。 定 位 结果 与 精 轨 坐标 差 的 比较 见 
图 12。 红 色 星 号 为 校正 前 结果 ， 绿 色 加 号 为 校正 后 结果 ， 蓝 色 圆 圈 为 新 处 理 程序 的 结果 。 经 过 基线 
和 台 站 时 延 校正 后 ， 卫 星 定位 精度 有 一 定 提高 ， 尤 其 是 对 图 6 中 标 出 的 3 处 存在 异常 值 的 时 间 段 。 用 
新 的 后 处 理 程序 重新 处 理 后 的 轨道 定位 结果 与 精 轨 比 较 的 均 方 根 误差 约 6.61 米 ，3 次 处 理 结果 的 残余 
时 延 残 差 均 值 和 卫星 坐标 均 方 根 误差 统计 见 表 2。 


Fig. 12 Comparison of Satellite Coordinate Differences 
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试验 结果 表明 ， 利 用 仅 4 小 时 的 短 时 VLBI 时 延 数据 可 以 得 到 外 符合 精度 米 级 的 GEO 卫星 定位 
结果 ， 有 利于 GEO 卫星 的 快速 轨道 恢复 。 本 次 实测 试验 结果 与 历史 实测 结果 对 比 见 表 3， 对 比 说 明 ， 
此 次 VLBI 罕 带 观测 时 延 测量 精度 较 历 史 罕 带 观测 结果 有 较 大 提高 ， 且 接近 宽带 观测 的 时 延 精度 结果 。 
并 且 本 文 首次 计算 了 外 符合 精度 结果 ， 更 好 地 评估 定位 结果 的 准确 度 。 虽 与 精 轨 结果 比 仍 有 一 定 差异 ， 
但 在 GEO 卫星 轨 控 后 可 快速 提供 较 高 精度 卫星 位 置信 息 ， 减 少 轨 道 误差 影响 ， 仍 具有 重要 意义 。 
表 3 GEO 卫星 实测 试验 信息 统计 
Tab. 3 Information Statistics of GEO Satellites Observations 


year VLBI observation mode VLBI observation accuracy estimated horizontal position accuracy external coincidence accuracy 
1984 17 hours, 4GHz 2ns / / 
3 hours, 2/8GHz, 4M 
2003 bandwidth ^ DOR 41cm 8m / 
3 hours, 2GHz, 2M 
2011 bandwidth 3.6ns 10m / 


36 hours, 4GHz, 32M i 
2017 bandwidth 0.14ns 0.42m / 


4 hours, 2GHz, M 


our article Ac aul 0.16ns 0.58m 6.61m 


GEO 卫星 的 短 时 定位 精度 ， 有 利于 卫星 机 动 后 的 轨道 快速 恢复 。 本 文 首 先 对 VLBI 时 


bandwidth 


利用 VLBI 技术 观测 GEO 卫星 ， 可 充分 发 挥 VLBI 对 卫星 横向 位 置 约束 精度 高 的 特点 ， 提 高 
nr 仑 精 


度 分 析 ， 包括 射电 源源 位 置 、 台 站 位 置 、 地 球 自转 参数 和 介质 等 带 来 的 系统 性 误差 ， 和 观测 源 噪 声 、 


氧 原子 钟 不 稳定 度 、 介 质 扰动 和 仪器 相位 抖动 等 带 来 的 随机 误差 。 分 析 得 差分 VLBI 时 
度 为 0.0720ns， 换 算 为 卫星 平 位 置 误差 约 0.26 米 。 


才 延 理论 测量 精 


实测 实验 利用 CVN 的 4 个 台 站 观测 GEO 卫星 ， 得 4 小 时 的 VLBI 残 余 时 延 数据 的 拟 后 残 差 中 误 


差 均 值 为 1.03ns。 画 图 结果 显示 残余 时 延 存在 缺 值 和 跳 变 的 情况 ， 影 响 了 时 延 测量 精度 。 因 此 本 文 提 
出 了 针对 该 情况 的 异常 值 检测 校正 方法 ， 分 别 考虑 与 观测 台 站 和 观测 基线 相关 的 异常 值 校 正 。 重 新 处 


理 校正 后 的 数据 得 残余 时 延 拟 后 残 差 中 误差 均值 为 0.16ns。 进 而 进行 卫星 定位 计算 进 一 步 验证 该 时 延 


异常 值 的 校正 效果 ， 发 现 卫星 定位 结果 与 IGS 提供 的 精 轨 结果 差异 由 16.76 米 提 高 到 6.61 米 ， 说 明了 


校正 方法 的 有 效 性 。 本 文 的 研究 可 利用 短 时 的 VLBI 观测 数据 进行 GEO 卫星 的 快速 高 精度 定位 ， 有 


EZ 


]T GEO 卫星 轨道 的 快速 恢复 。 
在 此 基础 上 ， 未 来 还 可 利用 VLBI 时 延 观测 数据 和 时 延 率 数据 预报 短 时 间 的 VLBI 时 延 值 ， 并 进 


行 后 续 GEO 卫星 位 置 预报 的 相关 研究 。 随 VLBI 全 球 观测 系统 的 不 断 完 善 ， 将 会 有 更 多 台 站 和 更 多 


基线 组 合 可 用 于 GEO 卫星 的 跟踪 观测 。 另 外 可 以 通过 设置 专门 的 DOR 信 标 或 宽带 信 标 等 方式 ， 进 


一 步 提高 短 时 测量 精度 ， 提 高 GEO 卫星 轨道 监测 能 力 。 
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Fast high accuracy positioning of GEO satellite in VLBI 
observation 
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Abstract: With the rapid development of satellites Earth observation, navigation and positioning, the 
orbit of Geostationary Earth Orbit (GEO) satellites requires an accuracy of meters level or even higher. In 
this paper, we study the method of GEO satellites fast positioning using Very Long Baseline 
Interferometry (VLBI) delay data. First of all, we analyze the theoretical accuracy of VLBI delay. After 
that we conduct the GEO satellite observation experiment and data correlation using Chinese VLBI 
Network. Finally, the positioning calculation of GEO satellite is completed by using short-time VLBI 
delay data. In order to solve the problem of outliers related to stations and baselines in VLBI delay, we 
also study the corresponding detection and correction methods. After correction and reprocessing, the 
measurement accuracy of VLBI residual delay is about 0.16ns, which is equivalent to the transverse 
position error of about 0.58 meters. The research in this paper can be used for fast and high-precision 
positioning and rapid orbit recovery of GEO satellites. 

Key words: differential VLBI, GEO satellite, accuracy analysis, delay correction, fast high accuracy 


positioning 
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